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Abstract — In this paper, we describe the spectral charac- 
teristics that can be achieved in fiber reflection (Bragg) and 
transmission gratings. Both principles for understanding and 
tools for designing tiber gratings are emphasized. Examples are 
given to illustrate the wide variety of optical properties that 
are possible in fiber gratings. The types of gratings considered 
include uniform, apodized, chirped, discrete phase-shifted, and 
superstructure gratings; short-period and long-period gratings; 
symmetric and tilted gratings; and cladding-mode and radiation- 
mode coupling gratings. 

Index Terms — Distributed feedback devices, gratings, optical 
fibers, optical fiber communication, optical tiber devices, optical 
fiber filters. 



I. Introduction 

THE fiber phase grating written by ultraviolet light into 
the core of an optical fiber has developed into a critical 
component for many applications in fiber-optic communication 
and sensor systems. These are described in detail in this issue, 
in numerous previous papers, and in several excellent review 
articles [l]-[6]. Advantages of fiber gratings over competing 
technologies include all-fiber geometry, low insertion loss, 
high return loss or extinction, and potentially low cost. But the 
most distinguishing feature of fiber gratings is the flexibility 
they offer for achieving desired spectral characteristics. Nu- 
merous physical parameters can be varied, including: induced 
index change, length, apodization, period chirp, fringe tilt, and 
whether the grating supports counterpropagating or coprop- 
agating coupling at a desired wavelength. By varying these 
parameters, gratings can be made with normalized bandwidths 
(AA/A) between 0.1 and 10""^, extremely sharp spectral 
features, and tailorable dispersive characteristics. 

In this paper, we focus on the optical properties of fiber 
phase gratings, in order to provide the reader with both intu- 
ition for understanding and tools for designing fiber gratings 
and related devices. Though many of the concepts developed 
below are dispersed elsewhere in the literature, it is hoped that 
by drawing these together and putting them specifically in the 
context of fiber gratings, the unfamiliar reader will benefit from 
a concise starting point while the experienced reader will gain 
a clearer understanding of optical characteristics achievable 
in fiber-grating components. The discussion below is limited 
to the linear optical properties of fiber gratings. Although a 
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Fig. 1. Effective index parameter h and core confinement factor F versus 
normalized frequency for the LPoi mode of a step index fiber. 



number of interesting nonlinear-optical properties have been 
reported recently (see [7] for an excellent review), these are 
outside the scope of this article. However, the description of 
linear properties of fiber gratings below is fundamental for 
understanding these nonlinear properties. 

Fiber phase gratings are produced by exposing an optical 
fiber to a spatially varying pattern of ultraviolet intensity. Here 
we assume for simplicity that what results is a perturbation to 
the effective refractive index ^efr of the guided mode(s) of 
interest, described by 



Crieff (^) = 8nQf[{z)i 1 + v cos 



27r 

T 



z + <t>{z) 



(1) 



where ^rieH is the ''dc" index change spatially averaged over 
a grating period, v is the fringe visibility of the index change, 
A is the nominal period, and (j){z) describes grating chirp. If 
the fiber has a step-index profile and an induced index change 
8nco{z) is created uniformly across the core, then we find 
that (Crieff ^ Vbuco where F is the core power confinement 
factor for the mode of interest. For example, Fig. 1 shows 
the confinement factor F and the effective index parameter b 
for the LPoi mode. For LP;^ modes, i is a solution to the 
dispersion relation [8] 



Ki {Vs/h) 



(2) 



where / is the azimuthal order of the mode and V = 
(27r/A)a\/n|~-^nJ^ is a normalized firequency, with a 
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Fig. 2. Common types of fiber gratings as classified by variation of the 
induced index change along the fiber axis, including (a) uniform with 
positive- only index change, (b) Gaussian-apodized, (c) raised-cosine-apodized 
with zero -do index change, (d) chirped, (e) discrete phase shift (of and 
(f) superstructure. 



the core radius, Tico the core index and tici the cladding 
index. The effective index is related to b through 6 = 

(^eff ~ ^ ^d)- ^^^^ y aii^ ^ known, T can 

be determined from 



Y2 



1 - 



Jf(V^VT^) 



. (3) 



The optical properties of a fiber grating are essentially 
determined by the variation of the induced index change 
5nefi along the fiber axis z. Fig. 2 illustrates some common 
variations that are discussed in this paper. The terminology 
used to describe these is given in the figure caption. For 
illustrative purposes the size of the grating period relative to 
the grating length has been greatly exaggerated. 

The remainder of this paper is organized as follows. In 
Section II, we examine the optical properties of uniform fiber 
gratings when coupling occurs between only two modes, 
since this case yields simple solutions for the reflection and 
transmission. The coupled-mode theory that is a vital tool for 
understanding gratings is introduced in this section, and some 
basic design rules-of-thumb applicable to all fiber gratings 
are established. In Section III, the analysis is extended to 
include two-mode coupling in nonuniform gratings, including 
apodized, chirped, phase-shifted, and superstmcture gratings. 
Numerical techniques for analyzing these are discussed, and 
examples are provided. In Section IV, we briefly consider the 
effects of grating tilt, or blazing, on the optical properties 
of Bragg gratings. Then, in Section V, we discuss the opti- 
cal properties that result from coupling to a multiplicity of 
cladding modes, and in Section VI coupling to the continuum 
of radiation modes is considered. Finally, concluding remarks 
are given in Section VII. 
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Fig. 3. The diffraction of a light wave by a grating. 

II. Two-Mode Coupling in Uniform Gratings 

A. Resonant Wavelength for Grating Diffraction 

In Section II, we focus on the simple case of coupling 
between two fiber modes by a uniform grating. Before we 
develop the quantitative analysis using coupled-mode theory, 
it is helpful to consider a qualitative picture of the basic 
interactions of interest. A fiber grating is simply an optical 
diffraction grating, and thus its effect upon a light wave 
incident on the grating at an angle Oi can be described by 
the familiar grating equation [9] 



n sin ^2 = ^ sin 9i + m 



X 
A 



(4) 



where Oo is the angle of the diffracted wave and the integer 
m determines the diffraction order (see Fig. 3). This equation 
predicts only the directions O2 into which constructive interfer- 
ence occurs, but it is nevertheless capable of determining the 
wavelength at which a fiber grating most efficiently couples 
light between two modes. 

Fiber gratings can be broadly classified into two types: 
Bragg gratings (also called reflection and short-period grat- 
ings), in which coupling occurs between modes traveling in 
opposite directions; and transmission gratings (also called 
long-period gratings), in which the coupling is between modes 
traveling in the same direction. Fig. 4(a) illustrates reflection 
by a Bragg grating of a mode with a bounce angle of Ox into the 
same mode traveling in the opposite direction with a bounce 
angle of O2 ~ -^i- Since the mode propagation constant p 
is simply = (27r/A)ue£E where n^ff = Uco sin 0, we may 
rewrite (4) for guided modes as 

/i2=/ii+m — . (5) 

For first-order diffraction, which usually dominates in a fiber 
grating, m = — 1. This condition is illustrated on the /3 axis 
shown below the fiber. The solid circles represent bound core 
modes {rici < n^^ < rico), the open circles represent cladding 
modes (1 < nefr < n^O, and the hatched regions represent 
the continuum of radiation modes. Negative values describe 
modes that propagate in the ~z direction. By using (5) and 
recognizing p2 < 0, we find that the resonant wavelength for 
reflection of a mode of index riej/^i into a mode of index 



(6) 
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Fig. 4, Ray-optic illustration of (a) core-mode Bragg reflection by a fiber 
Bragg grating and (b) cladding-mode coupling by a fiber transmission grating. 
The p axes below each diagram demonstrate the grating condition in (5) for 
rn ~ ~1. 



If the two modes are identical, we get the famiUar result for 
Bragg reflection: A = 2neffA. 

Diffraction by a transmission grating of a mode with a 
bounce angle of Oi into a co-propagating mode with a bounce 
angle of Q2 is illustrated in Fig. 4(b), In this illustration the 
first mode is a core mode while the second is a cladding mode. 
Since here 62 > 0, (5) predicts the resonant wavelength for a 
transmission grating as 



(7) 



For copropagating coupling at a given wavelength, evidently 
a much longer grating period A is required than for counter- 
propagating coupling. 

B. Coupled-Mode Theory 

Coupled-mode theory is a good tool for obtaining quantita- 
tive information about the diffraction efficiency and spectral 
dependence of fiber gratings. While other techniques are 
available, here we consider only coupled-mode theory since 
it is straightforward, it is intuitive, and it accurately models 



the optical properties of most fiber gratings of interest. We 
do not provide a derivation of coupled-mode theory, as this is 
detailed in numerous articles and texts [10], [11]. Our notation 
follows most closely that of Kogelnik [1 1]. In the ideal-mode 
approximation to coupled-mode theory, we assume that the 
transverse component of the electric field can be written as a 
superposition of the ideal modes labeled j (i.e., the modes in 
an ideal waveguide with no grating perturbation), such that 



' Cjt{x^ y) exp(— zojt) 



(8) 



where Aj(z) and Bj{z) are slowly varying amplitudes of the 
jth mode traveling in the +z and —z directions, respectively. 
The transverse mode fields ejt{x, y) might describe the bound- 
core or radiation LP modes, as given in [8], or they might 
describe cladding modes. While the modes are orthogonal in 
an ideal waveguide and hence, do not exchange energy, the 
presence of a dielectric perturbation causes the modes to be 
coupled such that the amplitudes Aj and Bj of the jth mode 
evolve along the z axis according to 



3 _ 



dz 



k 

+ i B^Ki^ - K^^^) exp + fij)z] 



(9) 



dz 



k 

- i ^ Bu{Kl^ + Kl^) exp - /i;)^]. 



(10) 



In (9) and (10), K\^(z) is the transverse coupling coefficient 
between modes j and k given by 



00 



(11) 



where As is the perturbation to the permittivity, approximately 
Ae ^ 2n8n when 6n n. The longitudinal coefficient 
K^j{z) is analogous to Klj{z), but generally K^j{z) < 
Kl.j{z) for fiber modes, and thus this coefficient is usually 
neglected. 

In most fiber gratings the induced index change Sn{x^ z) 
is approximately uniform across the core and nonexistent 
outside the core. We can thus describe the core index by 
an expression similar to (1), but with ^nefr(^) replaced by 
bucoiz). If we define two new coefficients 

^kji^^) - ^^^-^ ^nco(^) / / dxdyeutix, y) - y) (12) 



core 



V 



(13) 



where a is a "dc" (period-averaged) coupling coefficient and 
t^ is an "AC" coupling coefficient, then the general coupling 
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coefficient can be written 

^hi^^) = <='kj{z) + 2n.kj{z) cos 



27r , , , 
-J-Z + <p{z) 



(14) 



Equations (9)-(14) are the coupled-mode equations that we 
use to describe fiber-grating spectra below. 

C Bragg Gratings 

Near the wavelength for which reflection of a mode of 
amplitude A{z) into an identical counter-propagating mode of 
amplitude B{z) is the dominant interaction in a Bragg grating, 
(9) and (10) may be simplified by retaining only terms that 
involve the amplitudes of the particular mode, and then making 
the "synchronous approximation" [11], The latter amounts to 
neglecting terms on the right-hand sides of the differential 
equations that contain a rapidly oscillating z dependence, 
since these contribute little to the growth and decay of the 
amplitudes. The resulting equations can be written 

= iaR(z)+inSiz) (15) 



dz 
dS 

dz 



~iaS{z) — iti*R{z) 



(16) 



where the amplitudes R and 5 are R{z) = A{z) exp {i6z — 
(/>/2) and S{z) ~ B{z) exp {~i6z + (/>/2). In these equations 
K is the "AC*' coupling coefficient from (13) and a is a general 
"dc" self-coupling coefficient defined as 



(17) 



The detuning S, which is independent of z for all gratings, is 
defined to be 



= 0-0D 

1 

= 27rnefi I ^ ~ ^ 



(18) 



where = 2neff A is the "design wavelength" for Bragg 
scattering by an infinitesimally weak grating (Crieff — ^ 0) with 
a period A. Note when 5 = 0 we find A = 2neff A, the Bragg 
condition predicted by the qualitative grating picture above. 
The "dc" coupling coefficient a is defined in (12). Absorption 
loss in the grating can be described by a complex coefficient 
a, where the power loss coefficient is a = 2 Im (a). Light 
not reflected by the grating experiences a transmission loss 
of lOlogj^Q (e)Q; dB/cm. Finally, the derivative {l/2)d4)/dz 
describes possible chirp of die grating period, where (/?{z) is 
defined through (1) or (14). 

For a single-mode Bragg reflection grating, we find the 
following simple relations: 



27T -T^ 



(19) 



If the grating is uniform along z^ then <^neff is a constant and 
d(j)/dz — 0, and thus a, and a are constants. Thus, (15) and 
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Fig. 5. Reflection spectra versus normalized vi^avelength for Bragg reflection 
in uniform gratings with kL — 2 (dashed line) and kL = 8 (solid line). 



(16) are coupled first-order ordinary differential equations with 
constant coefficients, for which closed-form solutions can be 
found when appropriate boundary conditions are specified. The 
reflectivity of a uniform fiber grating of length L can be found 
by assuming a forward-going wave incident from z = —oo 
[say R{~L/2) ~ 1] and requiring that no backward-going 
wave exists for z > L/2 [i.e., S{L/2) = 0]. The amplitude 
and power reflection coefficients p = S{—L/2)/R{—L/2) and 
r = lycp, respectively, can then be shown to be [10], [11] 



P ^ 



—K sinh (V^^ — a^Lj 



a sinh [V — a^L) + iV/s:^ — a'^ cosh (V k.-^ — a^L) 

(21) 



and 



r — 



cosh' {Vk^ - aU) - 



9 

a" 



(22) 



A number of interesting features of fiber Bragg gratings 
can be seen from these results. Typical examples of the power 
reflectivity r for uniform gratings with nL = 2 and kL ~ 8 
are shown in Fig. 5, plotted versus the normalized wavelength 

' ' (23. 



A 



max 



1 + 



where is the total number of grating periods {N =^ L/A), 
here chosen to be N — 10000, and A^ax is the wavelength 
at which maximum reflectivity occurs. If were larger 
or smaller, the reflection bandwidth would be narrower or 
broader, respectively, for a given value of kL. From (22), 
we find the maximum reflectivity r'j^ax for a Bragg grating is 



- tanh^ {kL) 



(24) 



(20) and it occurs when <j = 0, or at the wavelength 



Amax = ( IH ~ ]^D- 

^eff 



(25) 
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The points on these plots denoted by open circles indicate the 
''band edges," or the points at the edges of the "band gap"; 
defined such that |a| < k. Inside the band gap, the amplitudes 
R(z) and S{z) grow and decay exponentially along z\ outside 
the band gap they evolve sinusoidally. The reflectivity at the 
band edges is 



'■bandedge - ^ _^ ^^j^y 

and the band edges occur at the wavelengths 

^bandedge — ^max ^ ^ ^P- 

2neff 



(26) 



(27) 



From (27) the normalized bandwidth of a Bragg grating as 
measured at the band edges is simply 



and edge 



A 



(28) 



where vdn^^ is simply the "AC" part of the induced index 
change. 

A more readily measurable bandwidth for a uniform Bragg 
grating is that between the first zeros on either side of the 
maximum reflectivity. Looking at the excursion of a from 
^ = 0 that causes the numerator in (22) to go to zero, we find 



AAo 
A 




(29) 



In the "weak-grating limit," for which vSn^^fi is very small, 
we find 



AA 



0 



Ad 



vdUe^fi < 



A 



D 



L 



(30) 



the bandwidth of weak gratings is said to be "length-limited." 
However, in the "strong-grating limit," we find 



AA, 



v6n 



eflf 



L 



(31) 



In strong gratings, the light does not penetrate the full length of 
the grating, and thus the bandwidth is independent of length 
and directly proportional to the induced index change. For 
strong gratings the bandwidth is similar whether measured at 
the band edges, at the first zeros, or as the full-width-at-half- 
maximum (FWHM). 

To compare the theory to an experimental measurement, 
Fig. 6 shows the measured (dots) and calculated (line) re- 
flection from a LO-mm-long uniform grating with a design 
wavelength of 1558 nm and an induced index change of 
v6n^fi = 8 X lO"'^, thus yielding kL ~ 1.64. 

Recently, there is growing interest in the dispersive proper- 
ties of fiber Bragg gratings for applications such as dispersion 
compensation, pulse shaping, and fiber and semiconductor 
laser components. Although many of these rely on the ability 
to tailor the dispersion in the nonuniform gratings discussed in 
Section III, here we introduce the basis for determining delay 
and dispersion from the known (complex) reflectivity of a 
Bragg grating. The group delay and dispersion of the reflected 
light can be determined from the phase of the amplitude 
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Fig. 6. Measured (dots) and calculated (line) reflection spectra for Bragg 
reflection in a l~mm-long uniform grating with kL = 1.64. 



reflection coefficient /? in (21). If we denote &p ~ phase{p), 
then at a local frequency too we may expand 0p in a Taylor 
series about ojq. Since the first derivative dOp/duis directly 
proportional to the frequency u, this quantity can be identified 
as a time delay. Thus, the delay time for light reflected off 
of a grating is 



de 



Tp^ 



A2 dO, 



(32) 

duj 2t\c a a 
Tp is usually given in units of picoseconds. Fig. 7 shows the 
delay calculated for the two example gratings from Fig. 5. 
Here the grating length is L — 1 cm, the Design wavelength is 
\d ~ 1^50 nm, neff — 1.45, and the fringe visibility is -u = 1. 
For the weaker grating in Fig. 7(a), bn^^ = 1 x 10"""^, while 
for the stronger grating in Fig. 7(b), hn^fi — 4 x 10"'*. We see 
that for unchirped uniform gratings both the reflectivity and 
the delay are symmetric about the wavelength Ai^ax- 

Since the dispersion dp (in ps/nm) is the rate of change of 
delay with wavelength, we find 



dp = 



dTp 
dX 

2Tn A^ d^e^ 



A 27rc dX^ 
27rcd^r. 



(33) 



A2 du;^ 

In a uniform grating, the dispersion is zero near Amax? and only 
becomes appreciable near the band edges and side lobes of the 
reflection spectrum, where it tends to vary rapidly with wave- 
length. Qualitatively, this behavior of the delay and dispersion 
(along with numerous other characteristics of fiber gratings) 
can be nicely explained by an "effective medium picture" 
developed by Sipe et al. [12]. For wavelengths outside the 
bandgap, the boundaries of the uniform grating (at z = ±L/2) 
act like abrupt interfaces, thus forming a Fabry-Perot-like 
cavity. The nulls in the reflection spectrum are analogous to 
Fabry-Perot resonances — at these frequencies light is trapped 
inside the cavity for many round-trips, thus experiencing 
enhanced delay. 
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Fig. 7. Calculated reflection spectra (dotted line) and group delay (solid line) for uniform Bragg gratings with (a) kL — 2 and (b) kL = S. 



Transmission Gratings 

Near the wavelength at which a mode "1" of amplitude 
Ai (z) is strongly coupled to a co-propagating mode "2" with 
amplitude A2{z), (9) and (10) may be simplified by retaining 
only terms that involve the amplitudes of these two modes, 
and then making the usual synchronous approximation. The 
resulting equations can be written 

= iaR(z) + inS{z) (34) 



dz 
dS 

dz 



= ~~id-S{z) -h in'R[z) 



(35) 



where the new amplitudes R and 5 are i^(^) = 
Aiex])[~i{aii + 022)^^ 1^ exp(i5;2 - (j)j2) and 
S{z) = A2exp[-i(aii + 0-22)^/2] exp(-i^;j + (j)/2), 
and where an and a 22 are "dc" coupling coefficients defined 
in (12). In (34) and (35), k = K21 is the "ac" 

cross-coupling coefficient from (13) and a is a general "dc" 
self-coupHng coefficient defined as 



<^i\ - 022 



1 d(j> 

2 Iz 



(36) 
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Fig- 8. Calculated cross transmission ty^ through uniform transmission gratings with kL = Tr/2 (dashed line) and kL = 37r/2 (solid line). 



Here the detuning, which is assumed to be constant along z, is 



— TfAn 



1 



eft 



A Ad 



(37) 



where A^ = ArteffA is the Design wavelength for an infinites- 
imal ly weak grating. As for Bragg gratings, <^ — 0 corresponds 
to the grating condition predicted by the qualitative picture of 
grating diffraction, or A — AneffA. 

For a uniform grating a and k are constants. UnUke for 
the case of Bragg reflection of a single mode, here the 
coupling coefficient n generally can not be written simply 
as in (20). For coupling between two different modes in 
both Bragg and transmission gratings the overlap integrals 
(12) and (13) must be evaluated numerically. (For untilted 
gratings closed-form expressions exist, but these still require 
numerical evaluation of Bessel functions.) Like the analogous 
Bragg- grating equations, (34) and (35) are coupled first-order 
ordinary differential equations with constant coefficients. Thus, 
closed-form solutions can be found when appropriate initial 
conditions are specified. The transmission can be found by 
assuming only one mode is incident fi"om z — — oo [say 
-R(O) = 1 and 5(0) — 0]. The power bar and cross trans- 
mission, = \R{zf/\R{Q)\' and tx - |5(^)|V^(0)|^ 
respectively, can be shown to be [10] 




COS" t V + a'^z 
1 



IH 



sm" 



1 



sm 



1 + 



a' 




(38) 



(39) 



Typical examples of the cross transmission for two 
uniform gratings of length L with kL = 7r/2 (dashed line), 
kL = 37r/2 (solid line), and an ~ <722 ~ 0 are shown in 
Fig. 8, plotted versus the normalized wavelength defined in 
(23). Here the total number of grating periods is = 100. If 
A" were larger or smaller, the bandwidth would be narrower or 
broader, respectively, for a given value of nL. The maximum 
cross transmission (which occurs when <3" = 0) is 



(40) 



and it occurs at the wavelength 



Am ax — 



1 - (ail - *^22) 



27r 



(41) 



For coupling between a core mode "1" and a cladding mode 
"2" [see Fig. 4(b)] with an induced index change in the core 
only, (711 = (7 from (19), and we find that (J22 < cni (since 
the core confinement factor for the cladding mode is small). 
Therefore, we can approximate (41) as 



Amax — 1 + 



efi 



An, 



(42) 



where we have assumed that ^neff , the induced change in the 
core-mode effective index, is much smaller than Aneif . Com- 
paring (42) to (25), we see that the wavelength of maximum 
coupling in a long-period cladding-mode coupler grating shifts 
(toward longer wavelengths) as the grating is being written 
Tieff/Aneff times more rapidly than the shift that occurs in a 
Bragg grating. 

A useful measure of the bandwidth of a transmission grating 
is the separation between the first zeros on either side of the 
spectral peak. Looking at the excursion of a from & = Q that 
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Fig. 9. Measured (dotted line) and calculated (solid line) bar transmission through a uniform cladding-mode transmission grating. 



causes the numerator in (39) to go to zero, we find 



2A 



1- 



kL 



TT 



(43) 



for a grating in which at most one complete exchange of power 
between the two modes occurs (kL < tt). For very weak 
gratings the normalized bandwidth is simply AAq/ A — 2/iV, 
just as for weak Bragg gratings. For strongly overcoupled 
gratings where kL > tt, the sidelobes become significantly 
more pronounced and hence, a better measure of the bandwidth 
is the FWHM of the envelope traced by the peaks of the 
sidelobes. Looking at the first factor in the expression for tx 
in (39), we find 



AA 



env 



2Xn 



A 



TrAn, 



eff 



(44) 



Since sidelobes are usually undesirable, most transmission 
gratings are designed such that k,L < 7r/2 where for the 
strongest gratings (large cross transmission) kL ^ 7t/2. 

To compare the theory to an experimental measurement, 
Fig. 9 shows the measured (dots) and calculated (line) bar 
transmission t~ for a relatively weak grating that couples 
the LP 01 core mode to the lowest-order cladding mode in 
a standard dispersion-shifted fiber (Anefi = 0,0042). The 
grating is 50 mm long and has a coupling-length product of 
kL = 0.39. 

Whereas, the bandwidth estimates above are generally ac- 
curate for practical Bragg gratings, the estimates (43) and 
(44) for transmission gratings can be poor if the effective- 
index dispersion is large near the resonant wavelength. The 
estimates assume Arieff is independent of wavelength- But 
if, for example, Afiefj oc A over an appreciable span of 
wavelengths, then from (37) the detuning 5 could remain small 
and thus allow strong coupling over the entire span, giving rise 



to a much broader grating! Dispersion is more of a concern 
for transmission gratings mainly because these tend to have 
broader bandwidths than Bragg gratings even in the absence 
of dispersion. 

While the results in this section apply rigorously only to 
coupling between two modes by a uniform grating, many of 
these "rules-of-thumb"; prove to be excellent approximations 
even for nonuniform gratings. 



IIL Two-Mode Coupling in Nonuniform Gratings 

In this section, we investigate the properties of nonuniform 
gratings in which the coupling occurs predominantly between 
two modes. We consider approaches to modeling both Bragg 
and transmission gratings and look at examples of several 
types of nonuniformity. 

Most fiber gratings designed for practical applications are 
not uniform gratings. Often the main reason for choosing 
a nonuniform design is to reduce the undesirable sidelobes 
prevalent in uniform-grating spectra; but there are many other 
reasons to adjust the optical properties of fiber gratings by 
tailoring the grating parameters along the fiber axis. It has been 
known for some time that apodizing the coupling strength of 
a waveguide grating can produce a reflection spectrum that 
more closely approximates the often-desired "top-hat" shape 
[13]"[15]. Sharp, well-defined filter shapes are rapidly becom- 
ing critical characteristics for passive components in dense 
wavelength-division multiplexed (DWDM) communications 
systems. Chirping the period of a grating enables the dispersive 
properties of the scattered light to be tailored [13]. Chirped 
fiber gratings are useful for dispersion compensation [16], for 
controlling and shaping short pulses in fiber lasers [17], and for 
creating stable continuous- wave (cw) and tunable mode-locked 
external-cavity semiconductor lasers [18], [19], Sometimes 
it is desirable to create discrete, localized phase shifts in 
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an otherwise periodic grating. Discrete phase shifts can be 
used to open an extremely narrow transmission resonance 
in a reflection grating or to tailor the passive fiher shape. 
Perhaps the most well-known application of discrete phase 
shifts is the use of a "quarter- wave" or tt phase shift in the 
center of a distributed-feedback laser to break the threshold 
condition degeneracy for the two lowest-order laser modes, 
thus favoring single-mode lasing [20], [21]. Recently interest 
has grown in gratings with periodic superstructure, in which 
the coupling strength, "dc" index change, or grating period 
is varied periodically with a period much larger than the 
nominal grating period A. Such ''sampled" gratings have been 
proposed and demonstrated for a number of applications [22], 
including use as a wavelength reference standard for DWDM 
systems [23]. Understanding effects of discrete phase shifts 
and superstructure has become critical recently with the advent 
of meter-long Bragg gratings for dispersion compensation 
produced by stitching together exposure regions formed with 
multiple phase masks [24], [25]. 

We consider two standard approaches for calculating the 
reflection and transmission spectra that result from two-mode 
coupling in nonuniform gratings. The first is direct numerical 
integration of the coupled-mode equations. This approach has 
several advantages, but it is rarely the fastest method. The 
second approach is a piecewise-uniform approach, in which 
the grating is divided into a number of uniform pieces. The 
closed-form solutions for each uniform piece are combined by 
multiplying matrixes associated with the pieces. This method 
is simple to implement, almost always sufficiently accurate, 
and generally the fastest. Other approaches are also possible, 
such as treating each grating half-period like a layer in a 
thin-film stack (Rouard's Method) [26]. Like the piecewise- 
uniform approach, this method amounts to multiplying a string 
of matrixes; but because the number of matrixes scales with 
the number of grating periods, this approach can become 
intractable for fiber gratings that are centimeters long with 
10^ periods or more. 

The direct-integration approach to solving the coupled-mode 
equations is straightforward. The equations have been given 
above: (15) and (16) apply to Bragg gratings and (34) and 
(35) are used for transmission gratings. Likewise the boundary 
conditions have been described above. For a Bragg grating of 
length L, one generally takes R{L/2) = 1 and S{L/2) = 0, 
and then integrates backward firom z = L/2 to z ~ —L/2, 
thus obtaining R{—L/2) and 5(— L/2). Since the transmission 
grating problem is an initial-value problem, the numerical 
integration is done in the forward direction from z — —L/2 
to z ~ L/2, starting with the initial conditions R(—L/2) — 1 
and 5 (—L/2) = 0, for example. Typically, adaptive-stepsize 
Runge-Kutta numerical integration works well. 

For modeling apodized gratings by direct numerical inte- 
gration, we simply use the ^-dependent quantities crf.j[z) and 
ill the coupled-mode equations, which give rise to a 
a{z) that also depends on z. For some apodized grating shapes, 
we need to truncate the apodization function. For example, 
fiber gratings are frequently written by a Gaussian laser beam, 



and thus have an approximately Gaussian profile of the form 

4 In 2z^ 



Suefiiz) = Sn^fi exp 



FWHM^ 



(45) 



where Sn^ff is the peak value of the "dc" effective index 
change and FWHM is the fiiU-width-at-half-maximum of the 
grating profile. Typically (45) is truncated at several times 
the FWHM, i.e., we choose L 3FWHM. Another common 
profile is the "raised-cosine" shape 



1 



1 + cos ( 

VFWHM/ 



(46) 



This profile is truncated at L = FWHM, where it is identically 
zero. Many other apodized profiles arc of interest as well, such 
as the "flat-top raised-cosine." 

Chirped gratings can be modeled using the direct integration 
technique by simply including a nonzero ;s:-dependent phase 
term {~)d(l)/dz in the self-coupling coefficient a defined 
in (17) and (36). In terms of more readily understandable 
parameters, the phase term for a linear chirp is 

V en L? ^^^^ 



2 dz 



dz 



where the "chirp" dX^/dz is a measure of the rate of change 
of the design wavelength with position in the grating, usually 
given in units of nanometers/centimeters. Linear chirp can also 
be specified in terms of a dimensionless "chirp parameter" F 
[13], given by 



F 



FWHM^ 



Z' 



~ — 47rn 



eff 



FWHM^ dXo 



A? 



D 



dz 



(48) 



is a measure of the fractional change in the grating period 
over the whole length of the grating. It is important to 
recognize that because chirp is simply incorporated into the 
coupled-mode equations as a ^-dependent term in the self- 
coupling coefficient a, its effect is identical to that of a 
"dc" index variation i7{z) with the same z dependence. This 
equivalence has been used to modify dispersion of gratings 
without actually varying the grating period [15], 

Incorporating discrete phase shifts and superstructure into 
the direct-integration approach is straightforward. For exam- 
ple, as the integration proceeds along z, each time a discrete 
phase shift is encountered a new constant phase shift is added 
in the expressions (1) or (14). In the coupled-mode equations 
(15) and (16), we thus multiply the current value of k by 
exp {i(p) where (p is the shifi; in grating phase. Superstructure 
is implemented through the z dependence in a{z) and k{z). 
For example, for sampled gratings we simply set = 0 in the 
nongrating regions. 

The often-preferred piecewise-uniform approach to model- 
ing nonuniform gratings is based on identifying 2x2 matrices 
for each uniform section of the grating, and then multiplying 
all of these together to obtain a single 2x2 matrix that 
describes the whole grating [27]. We divide the grating into 
M uniform sections and define Ri and Si to be the field 
amplitudes after traversing the section i. Thus for Bragg 



1286 



JOURNAL OF LIGHTWAVE TECHNOLOGY, VOL. 15, NO. 8, AUGUST 1997 



gratings we start with Rq = R{L/2) ~ 1 and 5o = S{L/2) = 

0 and calculate R{~L/2) = R^4 and S{-L/2) = 5m, while 
for transmission gratings we start with Ro — R{~L/2) = 1 
and Sq = 5(— L/2) — 0 and calculate R{L/2) = Rm and 
S{L/2) = Sm' The propagation through each uniform section 

1 is described by a matrix defined such that 



Ri 
Si 



(49) 



For Bragg gratings the matrix Ff is given by (50) shown at the 
bottom of the page where Az is the length of the ith uniform 
section, the coupling coefficients a and k are the local values 
in the ith section, and 



a' 



(51) 



Note 7b is imaginary at wavelengths for which \a\ > k. For 
transmission gratings the matrix FJ is shown in (52) at the 
bottom of the page where in this case 



7c 



(53) 



Once all of the matrices for the individual sections are known, 
we find the output amplitudes firom 



Ro 
So 



F — F^ • Fjv/— 1 * ■ . • * F^; ' . , . • Fi . 

(54) 



The number of sections needed for the piecewise-uniform 
calculation is determined by the required accuracy. For most 
apodized and chirped gratings M ~ 100 sections is sufficient. 
For quasiuniform gratings hke discrete-phase-shifted and sam- 
pled gratings, M is simply determined by the number of actual 
uniform sections in the grating. M may not be made arbitrarily 
large, since the coupled-mode-theory approximations that lead 
to (15), (16), (34), and (35) are not valid when a uniform 
grating section is only a few grating periods long [27]. Thus, 
we require Az ^ A, which means we must maintain 



(55) 



To implement the piecewise-uniform method for apodized 
and chirped gratings, we simply assign constant values a, ft, 
and {\j2)d(jy/dz to each uniform section, where these might be 
the ^-dependent values of <y{z), k{z), and {^)d4>/dz evaluated 



at the center of each section. For phase-shifted and sampled 
gratings, we insert a phase-shift matrix Fj^^ between the factors 
Fi and F^+i in the product in (54) for a phase shift after the 
ith section. For Bragg gratings the phase-shift matrix is of the 
form 



F^- = 



exp 



0 



0 



exp 




(56) 



— (F^J*, since we propa- 



and for transmission gratings F^^ — 
gate the field amplitudes in the ~z direction for Bragg gratings 
but the -^z direction for transmission gratings. Here, is the 
shift in the phase of the grating itself for discrete phase shifts, 
and for sampled gratings [see Fig. 2(f)] 

2 " A 



(57) 



where A^o is the separation between two grating sections. 

Having described two basic approaches for calculating 
reflection and transmission spectra through nonuniform grat- 
ings, we now give several examples that demonstrate the 
effects of apodization, chirp, phase shifts, and superstructure 
on the optical properties of fiber gratings. For most of the 
examples in the remainder of this section the piecewise- 
uniform method was used because of its speed, but the results 
are indistinguishable from those obtained by direct integration. 

To demonstrate the effects of apodization, Fig. 10 shows 
the reflection and group delay versus wavelength for gratings 
similar to those described in Fig. 7, only here the gratings have 
a Gaussian profile as illustrated in Fig. 2(b) and described 
by (45). The maximum index change values Suq^ and the 
FWHM — 1 cm for the gratings in Fig. 10(a) and (b) are 
the same as the uniform index change values and the length 
of the uniform gratings in Fig. 7(a) and (b), respectively. We 
see that the spectra are similar, except there are no sidelobes 
on the long-wavelength side and very different sidelobes on 
the short-wavelength side of the Gaussian spectra. The short- 
wavelength structure is caused by the nonuniform "dc" index 
change and has been described in detail elsewhere [12], [28]. 
Short wavelengths lie inside the local band gap (la| < k) 
associated with the wings of the grating and thus are strongly 
reflected there, but lie outside the local band gap near the 
center of the grating where they are only weakly reflected; the 
wings of the grating thus act like a Fabry-Perot cavity at short 
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Fig. 1 0. Reflection and group delay versus wavelength for Gaussian gratings 
similar to the uniform gratings in Fig. 7: (a) kFWHM = 2 and (b) 
/i-FWHM = 8. 
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wavelengths. Note that the difference in nonresonant delay (at 
wavelengths far from the Bragg resonance) between Figs. 7 
and 10 is not significant. The relative delay is more relevant 
than the absolute delay, which is sensitive to the identification 
of time zero (or ^ = L/2 in this case). 

In Fig. 1 1 (a) the reflection spectra for several nonuniform 
gratings are plotted, where for each the ''dc" index change 
is assumed to be zero (a 0), as illustrated in Fig. 2(c). 
In the models described above the "dc" index change is set 
to a small value (5nefr ^ 0) while the "ac" index change 
is maintained at the desired value, here vSn^^ = 1 x 10~^, 
by assuming a large fringe visibility v ^ oo. The grating 
profiles include a Gaussian, a raised-cosine, and a flat-top 
raised-cosine (in which the length of the uniform region 
is twice the FWHM of the cosine region) with the total 
FWHM — 10 mm for each. The Gaussian and raised-cosine 
gratings are indistinguishable on this plot. Note that the 
side lobes on the short-wavelength side of the Gaussian and 
raised-cosine spectra (see Fig. 10) have been eliminated by 
making the "dc" index change uniform, and the spectra are 
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Fig. 11. Calculated reflectivity of uniform and nonuniform gratings with 
"ac" index change of 1 X 10~^, zero "dc" index change, and FWHM =10 
mm. (a) Gaussian (dotted line), raised-cosine (solid line), and flat-top 
raised-cosine (dashed line), (b) uniform grating, and (c) expanded view 
of filter edge in dB. 
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Fig. 13. Calculated reflectivity of a uniform grating with periodic superstruc- 



Fig. 12. Calculated reflectivity of a raised-cosine grating with (solid line) (^^^^^ ^^^^) ^^^^^ 100-mm-long grating has 50 200-//in-long grating 

and without (dotted line) a tt phase shift in the phase of the grating at the sections spaced 1,8 mm apart. Also shown are the reflectivity of a single 
center. 200-/im-long section (dashed line) and of the compressed grating (dotted 

line), in which the 1.8 mm grating- free sections are removed. 



symmetric about A„aax. For comparison Fig. 11(b) shows a 
10-mm-long uniform grating with zero "dc" index change. 
To investigate how well the apodized gratings approximate 
a *'top-haf' spectrum. Fig. 1 1(c) shows a close-up view of the 
filter edges on a logarithmic scale. The (minimal) differences 
between the Gaussian and raised-cosine gratings are evident 
here; although the inherently truncated raised-cosine does 
exhibit side lobes, they are below —50 dB. 

An example of gratings with (solid line) and without (dotted 
line) a discrete phase shift are shown in Fig. 12. The gratings 
have a raised-cosine shape with FWHM = 10 mm, maximum 
"ac" index change vSn^fi = 2 x 10~^, and no "dc" index 
change (Crieff 0). The tt phase shift of the grating phase 
at the center opens a narrow transmission resonance at the 
design wavelength, but also broadens the reflection spectrum 
and diminishes r„iax- 

The reflection spectrum of a sampled grating exhibiting 
periodic superstructure is shown in Fig. 13 (solid line). As 
illustrated in Fig. 2(f), the grating is uniform except for regions 
where vSn^f^ is set to zero. Here, the "ac" index change in 
the nonzero regions is vSrieff ~ 5 x 10"'^, the total length 
of the grating is 100 mm, and there are 50 sections with an 
"on-off duty cycle of 10%. To facilitate an understanding 
of the resulting spectrum, the dashed line in Fig. 13 shows 
the reflection spectrum produced by a single section of this 
grating with the same index change; since its length is 200 
fj> m, kL ~ 0.2 for this section. The zeros in the envelope 
of the sampled-grating spectrum coincide with the zeros of 
the single-section spectrum. Also shown in this figure (dotted 
line) is a plot of the spectrum that results when the sampled 
grating is "compressed," such that the total length of the 
modulated portion of the grating is maintained (10 mm), but 
the unmodulated regions are removed. 

An example of a chirped grating is shown in Fig. 14. The 
grating has a "zero-dc" raised-cosine profile with FWHM ~ 1 
cm, vSn^fi = 5 X 10""^, and a chirp of (iXo/dz = — 1 
nm/cm. The plot shows the group delay calculated from 
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Fig. 14. Calculated group delay (dashed line) and dispersion (solid line) of 
a raised-cosine grating with "ac" index change of 5 x lO"'^, zero "dc" index 
change, FWHM =10 mm, and a chirp of — 1 nm/cm. The inset shows the 
reflectivity spectrum. 



(32) (dashed line), the dispersion calculated fi-om (33) (solid 
line), and the reflection spectrum (inset). The bandwidth of 
a similar, unchirped grating estimated from (31) is 0.53 
nm. A commonly used estimate of dispersion in a linearly 
chirped grating when variation of grating period along z is the 
dominant source of chirp is 



dp 



100 



-1 



D 



dz 



(ps/nm) 



(58) 



where dX^/dz is in imits nm/cm, and we have approximated 
2neff/c = 100 ps/cm. According to (58), the grating in Fig. 14 
should exhibit a dispersion of dp = —100 ps/nm and thus a 
delay of = 200 ps between 1549-1551 nm. These simple 
estimates disagree with the actual values by almost a factor 
of two, mainly because the effective chirp resulting fi^om 
apodization is comparable to the grating-period chirp. 



ERDOGAN: FIBER GRATING SPECTRA 



1289 




40- 



Uniform 
Gaussian 
Chirped Giuissi'rtn 



1549.5 



1550.0 

Wavelength (nin) 



1550.S 



155K0 



Fig. 15. Calculated bar transmission t= through a uniform grating (solid 
line), a Gaussian grating (dashed line), and a chirped Gaussian grating (dotted 
line) with grating parameters given in the text. 



Several examples of bar-transmission spectra through 
nonuniform transmission gratings are shown in Fig. 15. The 
solid line shows transmission through a uniform grating with 
a couphng-length product of kL ~ 7r/2, "zero-dc" index 
change (an ~ 022 = 0), an effective-index difference of 
Aneff = 0.1, and a length of 10 cm. The dashed line represents 
the transmission through a similar grating but with a Gaussian 
profile where the coupling-length product has been slightly 
reduced to /^-FWHM = 0.487r to maintain maximum coupling. 
Note the reduction of the side lobes> just as occurs for Bragg 
gratings. Although no example is shown, if the "dc" index 
change were nonzero, as shown in Fig. 2(b), the Gaussian 
grating spectrum would exhibit "Fabry-Perot" structure on the 
short-wavelength side of the spectrum similar to that shown 
in Fig. 10. The dotted line in Fig. 15 is the transmission 
spectrum obtained for a similar Gaussian grating with a linear 
chirp of d\u/dz = -0,025 nm/cm and • FWHM — 1.47r. 
Fig. 16 shows the spectrum of the group delay and dispersion 
associated with this grating; the cross transmission is shown 
in the inset. Although the dispersion is relatively small for this 
example, note how linear the delay is over the full bandwidth 
of the grating. 

IV. Tilted Gratings 

Another parameter the designer of fiber gratings has at his 
disposal is the tilt of the grating fringes in the core of the fiber 
[29], [30]. The main effect of grating tilt in a single-mode 
Bragg grating is to effectively reduce the fringe visibility v 
defined in (1), (13), and (14). However, grating tilt can also 
dramatically affect the coupling to radiation modes, which 
we discuss briefly in Section VI, as well as enable otherwise 
unallowed coupling between discrete bound modes of fiber. 
The latter interaction has been demonstrated in both Bragg 
[31] and transmission gratings [32] that couple the LPoi mode 
to the LP 11 mode of a dual-mode fiber. 

To see how grating tilt affects single-mode Bragg reflection 
in a fiber grating, suppose that the induced index change in the 
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Fig. 16. Group delay and dispersion for the cross amplitude of light trans- 
mitted through the chirped Gaussian grating described in Fig. 15. Cross 
transmission tx is shown in the inset. 
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Fig. 1 7. Diagram of the parameters associated with a tilted phase grating in 
the core of an optical fiber. 

core of the fiber 6nco is rotated by an angle 9, such that it is 

27r 



SUcoi'^^ ~ ^Tico{,z'~) \ 1 + COS 



V^9 



(59) 



where the z'-axis, shown in Fig. 17, is defined to be rJ ~ 
X sin 6 z cos 0. The grating period along the z (fiber) axis, 
which determines the resonant wavelengths for coupling, is 
thus Ag/ cos 9. For the slowly varying functions Snco{z^) 
and (}>{z^), we take z^ = z cos 0; i.e., we simply take the 
projection of these flinctions onto the fiber axis. The general 
coupling coefficient in (14) then becomes 

r27r 

^^iC-^) — ^i^) +2/v=p±(^) cos 



A 



z + <f}{z cos 0) 



(60) 



where we recognize that the subscripts j and k describe the 
same mode, except when one is associated with the forward- 
going mode (+) the other describes the backward-going mode 
(— ). The self and cross coupling coefficients now become 



cr(^) Sucoiz COS 0) j j dx dye^t{^, y) 




core 



and 



(61) 



Kzf±{z, 0) ~ -z — r — Snco{z COS 0) dxdy 



2 2 
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' exp ^±'i ^ X tan 0^ ^^(a;, y) 
'C*^^{x,y), (62) 
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Fig. 18. Plot of the normalized effective fringe visibility associated with 
single-mode Bragg reflection in a tilted grating in the core of a fiber with 
parameters described in the text. 

Notice that = {n^^y. Except for scaling of the slowly 
varying functions, the effects of tilt can be lumped into an 
"effective fringe visibility" v^±{d), defined such that 

V 

core 

j jdxdye^t{x, y) • ej.^ {x, y) 

core 

(63) 

Therefore, we may write 

n^±{z,0) = '^S^a[z) (64) 

in direct analogy to (13). This important result states that the 
effect of grating tilt on single-mode Bragg reflection is simply 
to reduce the effective fringe visibility by an amount given in 
(63). As an example, Fig. 18 shows a plot of the normalized 
effective fringe visibility as a function of tih angle 9 for the 
LPoi mode in a fiber with a cladding index of 1.44, a core- 
cladding A ~ 0.0055, a core radius of a — 2.625 /xm, and at 
a wavelength of 1550 nm. 

The Bragg reflection spectrum for a tilted grating can still 
be calculated from (15) and (16) where now k ~ 
Fig. 19 shows reflection spectra over a range of tilt angles 
for a Gaussian grating with a FWHM = 5 mm, a maximum 
index change of 8n^fi — 7,25 x 10"'^, and zero-tilt visibility 
V — 1. Because tilt affects Bragg reflection by effectively 
reducing the fringe visibility, this plot also demonstrates how 
the reflection spectrum of a Gaussian grating changes as the 
visibility is reduced (here by keeping the "dc" index change 
constant). Proper use of tilt for reducing the effective visibility 
and hence, Bragg reflection can be practically useful in cases 
where large radiation or other bound mode coupling is desired 
but Bragg reflection is undesirable [31]. 
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Fig, 19. Calculated reflectivity spectrum over a range of grating tilt angles 
for a Gaussian grating. 

V. Coupling to Cladding Modes 

There are numerous applications of fiber gratings that in- 
volve coupling of light into and out of the core mode(s) of 
a fiber. Fig. 20 shows measured spectra of the transmission 
of an LPoi core mode through an untilted Gaussian fiber 
grating that is 5 mm long and has a maximum index change 
of ^eff = 2 X 10"^. In Fig. 20(a) the bare (uncoated) 
section of fiber that contains the grating is immersed in 
index-matching fluid to simulate an infinite cladding. What 
results is a smooth transmission profile for A < 1540 nm 
demonstrating loss due to coupling of the core mode to the 
continuum of radiation modes. Radiation mode coupling is 
briefly discussed in Section VI and is described in detail 
in [28] and [29]. In Fig. 20(b), the fiber is immersed in 
glycerin, which has a refractive index greater than the cladding 
index. The transmission spectrum now exhibits fringes that are 
caused by Fabry-Perot-like interference resulting from partial 
reflection of the radiation modes off of the cladding-glycerin 
interface. In Fig. 20(c) the bare fiber is surrounded by air; 
for the wavelengths shown in this plot, the LPoi core mode 
couples into distinct backward-going cladding modes with a 
well-defined resonance peak for each mode. Here we briefly 
examine this last type of interaction. A more detailed descrip- 
tion can be found in [33]. 

Fig. 4(b) illustrates coupling of a core mode to a cladding 
mode from a ray point-of-view, where cladding modes are 
bound modes of the total fiber with 1 < nefi < Ud (assuming 
the fiber is surrounded by air). It is shown above that the 
grating period required for co-propagating coupling between a 
core mode of index n^^ and a cladding mode of index n^jj is 
A = ^/(^eff~''^efr)- analogous fashion, counter-propagating 
modes [^2 < 0 in Fig. 4(b)] can be coupled by a grating with 
a much shorter period A = A/(n^^ + rf^^). 

Using the coupled-mode theory outlined above, the equa- 
tions that describe coupling among the LPqi core mode and 
the (exact) cladding modes of order 1^, by a reflection grating 
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Fig. 20. Measured transmission through a Bragg grating where (a) the 
uncoated fiber is immersed in index-matching liquid to simulate an infinite 
cladding, (b) the fiber is imnnersed in glycerin, and (c) the bare fiber is 
surrounded by air and thus supports cladding modes. 



are [33] 

dz 



(65) 



E 



dB 



dz 



(66) 



where the detuning is 



= (67) 

Here, o-qIZ^i is the self-coupling coefficient for the LPqi 
mode given by (19), tt^i^^oi cross-coupling coefficient 

defined through (12) and (13), and we have neglected terms 
that involve self and cross coupling between cladding modes 
since the associated coupling coefficients are very small, or 

'^iM-t- /^t-L « « c^oili- In (65) and (66), we 

have also not included Bragg reflection of the LPoi core mode, 
making these equations valid only at wavelengths far away 
from this Bragg resonance. 
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Fig. 21. (a) Measured and (b) calculated transmission through a strong 
Gaussian Bragg grating, showing both core-mode Bragg reflection into 
cladding modes. 



The common synchronous approximation was used in ob- 
taining (65) and (66), but unlike (15) and (16) the resulting 
equations contain coefficients that oscillate rapidly along z. 
This is because there are now multiple values of S, which 
are not easily removed by defining R and S as we did for 
(15) and (16). However, if the cladding mode resonances do 
not overlap, then we may calculate each resonance separately^ 
retaining only the core mode and the appropriate cladding 
mode for the resonance. In that case, the problem reduces to 
simple two-mode coupling, as described in detail in Sections II 
and III. 

As an example of cladding-mode coupling in a Bragg 
grating, Fig. 21 shows the measured (a) and calculated (b) 
transmission spectrum through a Gaussian grating written in 
hydrogen-loaded [34] Coming Flexcore fiber. The grating 
length is FWHM — 4 mm, the index change is Suefi — 
2.8 X 10""^, and the visibility is approximately v ~ 1. The 
large transmission dip at 1541 nm is the Bragg reflection reso- 
nance^ while the remaining dips are cladding mode resonances. 
Here the resonances overlap substantially over most of the 
wavelength range shown in the plot. As a result, an accurate 
calculation could be obtained only by including multiple (but 
not necessarily all) cladding modes in (65) and (66) at each 
wavelength. The resolution in the experimental measurement 
is 0.1 nm. 

The equations that describe coupling among the LPqi core 
mode and the cladding modes of order \ i^hy a transmission 
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Fig. 22. Calculated transmission spectra through (a) a relatively weak and 
(b) a stronger transmission grating, each designed to couple the LPoi core 
mode to the jii = 7 cladding mode at 1550 nm. Solid lines represent a uniform 
grating while the dashed line represents a Gaussian grating. 



grating are [33] 



dz 



(68) 



E 



cl—co 



dz 



(69) 



Oi = - [ Poi — Piu, — 



27r 

A 



(70) 



In arriving at (68) and (69), approximations were applied 
similar to those made for (65) and (66). Fig. 22 shows calcu- 
lated examples of cladding-mode coupling in a transmission 
grating. Fig. 22(a) shows the transmission for a relatively weak 
grating with a length (FWHM) of 25 mm, a maximum index 
change of dn^fi = 2 x 10"^, and a uniform (solid line) or 
Gaussian (dashed line) profile. The four main dips seen in 
these spectra correspond to coupling to the = 1, 3, 5, 
7, and 9 cladding modes. Fig. 22(b) shows the transmission 
for a stronger, uniform grating with an induced index change 
of 1 X 10~^. In each case the grating period is adjusted to 
achieve coupling at 1550 nm between the LPqi core mode and 
the /i = o cladding mode. That is, the gratings in Fig. 22(a) 
and (b) have periods of A = 545 fjtm and A = 502 /rni, 
respectively. 



In the above discussion we consider only the Ip. cladding 
modes with azimuthal order / = 1. In an untilted grating 
coupling may occur only between modes of the same az- 
imuthal symmetry (as shared by the LPqi = HEn mode 
and the exact }{E/ERi^ cladding modes). In a typical 125- 
}jm diameter fiber there are several hundred ^ = 1 cladding 
modes at IR communications wavelengths. In a tilted grating, 
in which coupling to an enormous number of higher azimuthal- 
order cladding modes is allowed, the resonances are so densely 
spaced that they are generally not well resolvable. 

VI. Radiation-Mode Coupling 

Although numerous applications are possible for fiber grat- 
ings as waveguide-grating couplers that couple free-space 
beams into and out of bound fiber modes, surprisingly little 
work has been reported in this area. However, the effect of 
radiation-mode coupling as a loss mechanism on core-mode 
transmission has been studied [28], [29]. We do not describe 
this situation in detail, but merely introduce the reader to 
the main concepts in order to put radiation-mode coupling 
in context with the interactions discussed above. 

Fig. 23 illustrates the coupling of an LPqi core mode to the 
continuum of backward propagating radiation modes from a 
ray point of view. The coupling strength to these high-order 
(small-^2) modes is quite small unless the grating is suitably 
blazed. But as Fig. 20(a) illustrates, if the fiber is immersed in 
an index matching fluid or recoated with a high-index polymer, 
bound cladding modes cease to exist and substantial coupling 
to low-order (|/?2l > 27r/A) radiation modes can occur near the 
Bragg reflection resonance even in untilted gratings. A useful 
quantity to be aw^are of is the wavelength Acut at which true 
radiation-mode (not cladding-mode) coupling "cuts" on and/or 
off. Assuming the core confinement factor for the radiation 
modes is much smaller than that of the bound mode of interest, 
then we find 



A 



cut 




n 



1 + 



(71) 



where n = rici for the case of an infinitely clad fiber with index 
rtci [see Fig. 20(a)], orn — 1 when the fiber is surrounded by 
air such that bound cladding modes may propagate. The 
sign in the first factor applies to reflection, while the " sign 
corresponds to forward coupling. The second factor describes 
the shift of A^ut with increasing "dc" index change. 

Consider the coupling of an LPoi core mode to backward- 
propagating radiation modes labeled LP-pp, where the discrete 
index p identifies the polariza tion and azimuthal order, while 

the continuous label p = y(27r/A)2n^; - /i^^ denotes the 
transverse wavenumber of the radiation mode (0pp is the 
usual axial propagation constant). Here we assume the fiber 
has a cladding of index rid and with an infinite radius. The 
coupled-mode equations for this case are then 

dA 



dz 



P 



PP 



dB 



VP _ 



dz 



(72) 
(73) 
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Fig. 23. Ray -optic illustration of core-mode coupling to a 
backward- traveling radiation mode by a Bragg grating. The B axes 
below the diagram demonstrates the grating condition in (5) for m — 1. 



where 



(74) 



and where A is the amplitude of the core naode, Bp^ are the 
ampHtudes of the (continuous spectrum of) radiation inodes, 
and the usual summation now includes an integral in (72). 
Also, o'oiloi is the LPoi self-coupling coefficient given by 
(19)j and ^oi~^^ - i^p^oiT is the cross-coupling coefficient 
defined through (12) and (13). By applying essentially a "first 
Bom approximation" to the core-mode amplitude A, we can 
obtain an approximate expression for Bpp firom (73). After 
inserting this result into (72) and performing the integral over 
p, it can be shown [28]^ [29] that the core mode amplimde 
approximately obeys an equation of the form 



dA 
dz 



.CO—CO A 

Ji — 



'^^01-01- 



E 



L P 



p 



PP l^rtt— CO [2 



A (75) 



where the term in square brackets is evaluated at = 
27r/A — Poi. Since this term is real, clearly it gives rise to 
exponential loss in the amplitude of the core mode. Notice 
the loss coefficient is proportional to the square of the cross- 
coupling coefficient, and hence, to the square of the induced 
index change. 

Although the analysis given here is simplified, it is possible 
to predict the radiation-mode coupling loss spectrum even 
when the grating is tilted following a similar development 
to that in Section IV. As might be expected, as the tilt 
angle is increased (the grating is blazed), light can be cou- 
pled more efficiently at smaller angles |^2|- But increasingly 
many radiation-mode azimuthal orders must be included in 
the summation in (75) to accurately model radiation modes 
propagating more normal to the fiber axis. As an example of 
some typical transmission spectra, Fig, 24 shows a calculation 
of the loss in transmission through the same Gaussian grating 
described in Fig. 19, only here Sn^ff = 1.5 x 10~^. Spectra 
are shown for a range of tilt angles between 0 -45^, where 
the grating period A along the fiber axis is kept constant. The 
LpQi mode incident on the grating is assumed to be polarized 
perpendicular to both the x and z axes (see Fig. 17). The 
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Fig. 24. Calculated loss in transmission through a Gaussian grating in a fiber 
with an infinite cladding over a range of grating tilt angles. 



peak at the longest wavelength results from Bragg reflection, 
whereas the loss at other wavelengths is due to radiation-mode 
coupling. As expected, the efficiency for coupling to smaller 
and smaller angles B2 (which occurs at shorter wavelengths) 
improves as the grating is increasingly blazed. 

VIL Conclusions 

We have presented principles and tools for understanding the 
spectra exhibited by fiber Bragg and transmission gratings. It 
is hoped that the reader less familiar with these exciting com- 
ponents now has a base with which to understand the current 
and potential applications of fiber gratings, and that the more 
experienced reader has benefited fi*om a concise description of 
multiple fiber-grating phenomena with a uniform notation. 

An important principle to take away from Section II is that 
the simple closed-form expressions for optical properties of 
uniform fiber gratings are fi-equently excellent approximations 
for nonuniform gratings as well. This section also uncovers 
some of the basic similarities and differences between Bragg 
gratings and transmission gratings. One important difference 
is that as a Bragg grating is made increasingly stronger (by 
enlarging Sn^^% the maximum reflectivity and bandwidth 
simply increase; in contrast, to obtain a strong resonance in a 
transmission grating requires hitting a "target" Sn^f^ (usually 
nL = 7r/2) to avoid over-coupling. 

By examining nonuniform gratings in Section III, we found 
that it is critical to take into account both apodization and 
chirp of the grating period when the dispersive characteristics 
of fiber gratings are important. We also saw that it is almost 
always advantageous to minimize the "dc" apodization, unless 
it is intentionally utilized for tailoring dispersion, for example. 
In Section IV, we found that tilted gratings enable coupling 
between modes that have dissimilar azimuthal symmetry, and 
that the effect of tilt on single-mode Bragg reflection can be 
regarded simply as a reduction in the effective fringe visibility 
of the phase grating. Finally, in Sections V and VI, we briefly 
explored coupling from a core mode to cladding and radiation 
modes. Although long-period cladding-mode coupling gratings 
have recently seen rapid development [35], there appears to 
be relatively little investigation into the apparently abundant 
possibilities for applications based on exchange of light in the 
core of a single-mode fiber with the readily accessible cladding 
and radiation modes. 
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